DESY-11-074 



O 
(N 



in 
o 



Recombination within multi-chain 
contributions in pp scattering 



J. B artels" and M.G.Ryskin^ 



^ ! "II. Institut fiir Theoretische Physik, Universitat Hamburg, Luruper 

Chaussee 149, D-22761 Hamburg, Germany 



^ Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg, 188300, 



Oh. 

r^, ! Russia 

QQ I Abstract 

I We investigate the evolution of multiple parton chains in proton-proton scattering and 



show that interactions between different chains may become quite important. 



1 Introduction 



In the last years it has become clear that multiple parton interactions play an important 
^ I role in hadron-hadron collisions at high energies [Ill2l[3llH|5lini[ZllH]. As a first step 
^ ■ these chains are modelled as a collection of single noninteracting chains (Fig.l). Each 
- - - chain follows the usual partonic DGLAP evolution, i.e. the ladders are in color singlet 
states, and the momentum transfer across the ladder is set equal to zero. Disregarding 
final state radiation and working in leading order only, the final state produced by k such 
chains consists of N = ni + n2 + ... + nk partons, and the cross section, da ~ \T2^N\'^dfl]sf, 
is described by a sum of squares, without any interference terms. 

The theory of multiparton (higher twist) evolution has been outlined in |9j. To leading 
order, the evolution is described by the sum over the pairwise interactions between two 
i-channel partons, and the evolution kernels are given by the nonforward DGLAP split- 
ting functions. Of particular importance is the small-x region where powers of In 1/x may 
compensate (and even overcome) the higher twist suppression. In this region, the domi- 
nant contributions are given by gluon ladders, and their evolution in Inl/x is described 
by the BKP equations [10] . In leading order this evolution is, again, described by the sum 
over parwise interactions between t-channel gluons. The evolution kernels are given by 
the nonforward BFKL-kernels. At small x, the leading logarithmic approximations of the 
two approaches - higher twist evolution in momentum scale or small-x evolution in In 1/x 
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Figure 1: Two noninteracting chains: (a) the cross section (energy discontinuity of 
the scattering amphtude T2^2), (b) redrawn as a (uncut) two-ladder exchange diagram. 
Within each chain, the boxes mark the hard subprocess with largest transverse momenta 
(production of dijets). 



- coincide in the so-called 'double logarithmic approximation' which samples powers of 

{as In l/xlnp^). 

For both evolution schemes the t-channel multiparton state is in a color singlet state 
and, as far as the total cross section is concerned, the total momentum transfer is set 
equal to zero. However, any subsystem consisting of two ^-channel gluons, in general, 
will have nonzero color quantum number and nonzero momentum transfer. Therefore, 
describing the evolution of a t-channel state consisting of, say, 2n gluons as the evolu- 
tion of n noninteracting color singlet ladders with zero momentum transfer represents an 
approximation whose validity deserves further investigation. 

In this note we study, as a first correction beyond the approximation of noninteracting 
ladders, a particular type of 'interactions' between two ladders which we illustrate in Fig. 2: 
Starting at the proton at the bottom of Fig.2, we first have two noninteracting color singlet 
ladders (denoted by the pairs of t-channel gluons (14) and (23)). At rapidity Y we allow 
for a 'recombination' of t-channel gluons: from now on we have the two color singlet pairs 
(13) and (24). In the following we will denote his transition by 'recombination vertex'. It 
introduces a correlation between the two ladders. However, it is important to note that, 
in the double log approximation, this kind of interaction between the two ladders still 
belongs to the leading logarithmic approximation: for each momentum integral we have 
a factor (a^ Inp^ In l/x). It is this kind interaction between the two ladders which we 
will study in the following, staying within the double logarithmic approximation of gluon 
ladders. Particular attention will be given to the possibility that this recombination of 
two chains takes place in the perturbative region, i.e. in the region of large transverse 
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Figure 2: Two recombinations within the two chains. 



momenta. 

Recently an important potential application of such recombination effects has been 
suggested. In their attempt to explain the ridge effect reported by the CMS group at the 
LHC [H] , it has been suggested p!2] that the observed long range rapidity correlation and 
azimuthal correlation can be explained, within the Color Glass Condensate framework, 
by a two-chain recombination which will be discussed in this paper. 

2 Two noninteracting ladders 

We begin with the double logarithmic approximation of the two-chain configurations 
shown in Figs.l and 2: we search for regions of integration where each closed momen- 
tum loop gives two logarithms, one in the transverse momentum and one in rapidity. 
We restrict ourselves to gluon ladders which, at small x, are known to give the largest 
contributions. We parametrize our momenta as 

k = xpA + yPB + k (1) 

where pa, Pb are the large momenta of the incoming protons A and B, resp., the momen- 
tum fractions x, y range between and 1, and k denotes the two dimensional transverse 
momentum. In the double logarithmic approximation, the BFKL kernel and the splitting 
function Pgg lead to the same answer (Fig.3): we can either start from the small-x limit 
which is described by the BFKL equation and then take the limit of strongly different 
momentum scales; alternatively, we can begin with the coUinear limit where the DGLAP 
equations apply and then take the limit of small x. For our purposes it will turn out that 
the approach based upon the BFKL equation is more suitable: it is the region where the 
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Figure 3: A single cell inside a ladder, below the produced pair of jets. 



logarithms in 1/x are slightly larger than the transverse momentum logarithms where the 
recombination effects become important. 

In the region of strongly ordered transverse momenta the color singlet BFKL kernel 
is approximated by 



k\q ~ kf + k'\q ~ k 
(k - k'Y 



l\2 



2g'k'{k' - q) 



(2) 



This approximation will be used in the following. 

We begin with the two noninteracting ladders shown in Fig.l. Our main focus is on 
the integration over the loop momentum q, and we consider a single cell with momentum 
k' inside one of the ladders below the produced pairs of jets. This cell is illustrated in 
Fig.3. Using (jj]) for the upper rung (and the corresponding expression for the lower rung), 
together with the gluon propagators for the vertical lines, one sees that the integration over 
the transverse momentum k' is logarithmic only if k' 3> q, i.e. q defines the momentum 
scale Qo where the k^ evolution along the ladders starts. On the other hand, the range 
of the integration over the momentum transfer along the ladders is determined by the 
size of the interaction region and by the correlation length of the initial gluons of the two 
ladders inside the proton: we denote this effective radius by R, and put Ql = As 
a result, the cross section for the production of two pairs of gluon jets (cf. Fig.l) will be 
of the form: 



da 



1 



1 



f{xi,pl)f{yi,pl)f{x2, pl)f{y2, pI) 



(3) 



where f{x,p'^) = xg{x,p^) denotes the gluon density with initial momentum scale Qg, the 
factor results from the integration over the loop momentum q, and the momentum 

factors 1 / {p\p1)'^ represent the two production vertices of the pairs of gluon jets (evaluated 
in the double logarithmic approximation). 

As an important feature of ([3]) we mention that, as long as interactions between 
different parton chains are not taken into account, the momentum transfer along a ladder 
is of the order of the initial momentum scale of the (multi)parton distribution. 
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Figure 4: Kinematics of a recombination (14)(23) — )■ (13)(24). 



3 Recombinations in a two-ladder configuration 

We now turn to the main topic of this paper, a study of the two recombinations shown in 
Fig. 2. We begin with the recombination vertex below the produced jets; the kinematics 
are illustrated in Fig.4,. Above the recombination, the two ladders formed by the lines 
(13) and (24) are in the color singlet configuration; below the color singlet ladders are 
formed out of (14) and (23). As a result, each color recombination is accompanied by the 
color suppression factor 

In order to have, for the momentum loops below the lower two rungs (connecting (14) 
and (23)), logarithmic momentum integrals, we need 

\K\AK\^\i\- (5) 

Similarly, in order to find logarithmic momentum integrals above the upper two rungs 
(connecting (13) and (24)), we need 

\ki\, 1^2! > \q\, (6) 

Finally, inside the k loop of Fig.4. we need 

1^11,1^2! > > 1^2!, |fc| > |Z|, |fci|, 1^2! > 1^21- (7) 

Using, for the upper two rungs, the approximations following from ([2]), and combining 
them with the propagators for the vertical gluon lines, we obtain: 



{2N,g 



2\2 



d'^k f d^q 2k{k-q) 2k{k - q) 



{2Tif J i2TT)^ k\q - ky k\q - k)^ 



In order to obtain a logarithmic integrals of k, we identify two regions of k and q: either 
\q\ ^ |fc or |q| ^ \q — k\. For each of the two cases, after averaging over the azimuthal 
angle, we arrive at the integrals 



2q^ I r (9) 

^ ' (27r)3g2 ] (2vr)3fc2 
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which gives the desired logarithmic integral in k (or [q — k)). 

For the two ladders below the recombination we derive, from the condition ([7]) and 
from the second integral in (Q, that the upper cutoff is given by g^. The lower cutoff is 
obtained by applying the discussion of section 2.1: the loop momentum I appears only 
inside the initial condition of the lower proton B, and it is restricted by the effective scale 
Qo The condition implies that the momentum q also defines the lower momentum 
cutoff for the ladders above the recombination vertex. 

In order to find the full dependence on g^, we need to consider the full diagram in 
Fig.2 0. The recombination vertex above the produced pairs of jets is analysed in the 
same way as the lower one. This leads to an expression similar to (j9]), i.e. the complete 
dependence in q is of the form: 

(10) 

The integral in q is dominated by small values. Since the momentum q defines the upper 
momentum cutoff, both for the two ladders below the lower recombination vertex and 
for the two ladders above the upper recombination vertex, we conclude that the infrared 
divergence of the q-integration destroys the ladders above and below the recombination 
vertices. The recombination vertices, therefore, are absorbed into the nonperturbative 
initial conditions. As far as the perturbative part is concerned, we are back to the two 
noninteracting chains of section 2.1. 

The situation changes if logarithms in rapidity become more important than those 
in transverse momentum, i.e. within the BFKL approach we move towards small x. In 
Fig.2, we replace the rungs by BFKL Green's functions. Re-drawing the diagram in an 
more suitable way, we arrive at Fig. 5: Let us first reformulate the result which we have 
just obtained in the double logarithmic approxmation. We have shown that, in order 
to find inside the BFKL Green's function the maximal number of transverse logarithms, 
the momentm transfer across the Green's function has to be smaller than the transverse 
momenta of the two gluons entering the Green' function at the low-momentum side. In 
Fig. 5. this says that I and I' have to be small, while for the g-loop we found the integral 
J dq^ /q^ which favors small values, too. 

In order to see how the appearance of large rapidity intervals changes this situation, 
let us use the following integral representation for the forward BFKL-Green's function: 



^ ' ' ' J 2Tii J 2TTi \k'^ J 



^ - x(/^,o) 



(11) 



where /i = ii^ + |, the integration contours in /i and u run parallel to the imaginay axis, 
n) is the BFKL eigenvalue function, and we have averaged over the azimuthal angles 
of fc, k' . Furthermore, we have kept only the leading term of the conformal spin, n = 0. 
From this representation one easily deduces the dependence upon Y — Y' and In ^: after 
the integration over u the saddle point analysis of the remaining /i -integral shows that, for 
large In the dominant contributions come from /i = 1 ^ and oj = 0{1), whereas 
for large Y — Y', one finds yU ~ | and small u = ujbfkl = ^^^s^^^. This observation 



^Otherwise the momentum q will run through the blob corresponding to the proton initial conditions 
and, like the momentum Z, it will be restricted by a low scale Qo- 
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Figure 5: Another way of drawing Fig. 2b 



has important consequences. Let us denote the rapidities of the two recombinations by 
Y' and F, resp.. Beginning with the BFKL amphtudes near the produced jets, large 
rapidity intervals Y' — Yi, Yi — Y, Y' — Y2, Y2 — Y increase the corresponding anomalous 
dimensions and hence favor larger values of the central loop momentum, q. To see this 
im more detail, we insert the integral representations for all the BFKL Green' function 
in Fig. 5. We arrive at 



dcr 



1 1 1 



1 



dYidY2d^p,d^P2 Rl Rl {p\Y {pif 



dfi' 
2Txi 



dfi 
2Td 



dfi'i f dfii f dfi2 f dfi 



2T[i / 27ri / 27rz / 27rz 



Pi 



Mi 



,(y'-yi)xM) 



dY / dY' 



Pi 



Ml 



,{yi-y)x(m) 



SI 

Ql 



SI 

Ql 



Pi 



M2 



o{y'-y2)x(p'2) 



pI 



M2 



,(^2-y)x(M2) 



(12) 



Here we have introduced another length scale, Rl: it results from the additional integrals 
(as compared to eq.(l3])) d'^l and rf^Z', which are restricted by the proton radius and by the 
correlation between the two chains inside the proton. As long as all BFKL amplitudes are 
in DGLAP regime, i.e. they are dominated by the logarithms in the transverse momenta, 
all /i variables are small, and we are in the situation which we have described above: the 
q integral is dominated by small values, and the transverse momentum logarithms inside 
those four BFKL amplitudes which are close to the protons are destroyed. If, however, 
the rapidity intervals become large and the BFKL amplitudes are in the small-x region. 
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yU values are close to |. As a result, in f lT^ the overall power of may increase and the 
dominance of the small-q^ region disappears. One easily sees that large rapidity intervals 
near the protons, Ytot — Y' and Y , tend to make /x and n' large and thus help to increase 
the overall power of q^. 

Let us see in more detail how this balance works. Defining in fll2p the phase function 

/i'l, /ii, /i2, /i2, yu) = 2 {{Ytot - y')x{^^') + In ^ + 



+ {{Y' - Y,)xifi[) + fi[ In ^) + ((Fi - F)x(Mi) + /ii In ^ 
+ ({Y' - Y,)xifx',) + /i'l In 41 + f (>^2 - Y)xif^2) + /ii In 4 



Of ' VV -^ //V \r ^ / ' r _L n 

+2fyx(/^)+/iln^), (13) 



the saddle points are determined from the conditions: 



= | = 2,V.)r + 2l„J. (14) 



and 



which lead to 



and 



= ^ = xVM)(V'i-l^) + ln^, (15) 



In^ 

X'(/^.) = (16) 



In^ 



x'(/ii,.) = -^rfy (17) 

Similar equations are obtained for the other fi variables. 

For a systematic analysis one first determines, for fixed Y, Y' and g^, the stationary 
points of the n variables and then finds the dominant values of the rapidities Y , Y' and 
of the momentum scale Ing^. As we have said before, if in ( IT6l) the evolution in rapidity 
dominates over that in momentum scale, the rhs becomes small. Since x'{^^) vanishes at 
= we have 



1 1 



In^ 







^^"2 x"{\)' Y ■ 

On the other hand, if in f|T7|) the interval in momentum evolution is larger than in rapidity, 
the rhs is large when /i is close to zero: 



Y — Y 
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Figure 6: Two different paths of evolution: (a) Normal path which disfavors recom- 
bination, (b) a path which supports a recombination. For simplicity, we have taken 
pI = pI = and Yi = Y2. 



with a = and x"il) = a28C(3) (here C(3) ^ 1.202 denotes the Riemann zeta func- 
tion) . 

The results can be illustrated in terms of evolution paths in Fig. 6. There is an infinite 
number of paths in the In fc^-?/ plane which connect the protons with the produced jets 
with rapidity Yi, Y2 and momenta pf, p^- The saddle point analysis determines the most 
probable path. Two examples are shown in Fig.6. Let us now consider a few special cases. 
For simplicity, we start with the symmetric choice Yi = Y2 and pi = = p^. In order to 
get a large we search for the situation where [!> Hi and h' > iJi[. We insert the saddle 
point values (fTSl) . (fT9l) into (fT3l) . and first look for the extrema with respect to Y and F', 
that is for the saddle point of the expression 



which belongs to the lower part of Fig. 5. To get (!20|) we have used the value of fig in ( ITSj) 
and neglected a weak dependence of /i^ ~ 1/2 on y coming from (1181) . Thus the typical 
value of Y is 

ln^(»Vo^) 

leading to 

Aiower ~ exp (^-^ [2 \n'/\p'/q') - \n\p'/q')] + 2x(/is)l^i) • (22) 

An analogous expression is obtained for the upper part of the amplitude in Fig. 5, A^p. 
Assuming that the total available rapidity interval Ytot and the sub-rapidities Yi, Y2 are so 
large that the saddle point position f|T8l) is close to 1/2, i.e. 28C(3)aYi >> \n{q'^/Ql). 
We put fig = 0.5 in (I2U1) . fl2^ . and we see that in f[T^ the integral over q'^ takes the form 

/ ^e^P (i^ [21n^/^(pVg^) - ^ApVq')]) ■ (23) 
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Here we have used the LO BFKL ratio a/xifJ-s) = ct/xi^-^) = 1/4 In 2. The integral 
has its saddle point at 

q'^ r^p^exp{-z^) (24) 

with 2 = 3/4+ v/9/16 + 0.51n2 ~ 1.7; exp(— 2^) ~ 0.055. The corresponding value 
of fii^s — -2/4 In 2 ~ 0.6 is not small enough to justify the approximate estimate (IT^ . 
We therefore conclude that the competition between the Y and dependence leads to 
a rapidity saddle point fl2T]) somewhere inside the available interval {Yi,0), which in its 
turn leads to a rather large fii^s violating the initial inequality /x, < fi. 

In general we can say that the ordering /Xj < /i^ leads to the opposite ordering of 
the intercepts, xif^s) < x{f^i,s)- Therefore, in ( |T2l) the dominant contribution to the Y 
integral comes from the region of small Y where the recombination vertex is close to the 
initial proton and far from the production vertex of the dijets. For a small Y value the 
anomalous dimension /i cannot be large, and the essential values are small as well. This 
means that we are back to the situation of non-interacting ladders, illustrated in Fig. 6a. 

Next we turn to the opposite case fii > fi which corresponds to xIa^j) < Now 
the dominant F- value is large and close to the rapidity of high Et dijets Yi. However, 
now the whole anomalous dimension in the behaviour 2(/i + fi') — jii— fi'i — /i2 — < 
is negative, and the g^ integral is dominated by a low g^-value. 

The most interesting possibility is to put the recombination vertices just as close 
as possible to the high Et dijet production matrix elements. In this case there is no 
BFKL or DGLAP evolution in the intervals between the produced pairs of jets and the 
recombination vertices. That is, in the centre of Fig.5, we simply delete the four 'BFKL 
blobs' nearest to the produced jet pairs. Correspondingly, in (fT2|) we eliminate the third 
line, togather with the integrations over Hi, fi2, l^i, At2- The rapidities Y, Y' are close to YJ, 
and the g^ integral takes the form 



dlnq 



2 



g4(M.+M^) ^ (25) 



g- 

where the saddle point values, fig and fi'^, follow from the condition (HM : 

2 



and 



= x'{fis)Y + \n^ (26) 



= xVJmo*-r) + ln^. (27) 



Their values are taken from f|T8l) : 



2 x"{i) y 

i.e. the integral over g^ receives its main contribution from g^ close to mm{pl, p^} This 
situation belongs to the evolution path shown in Fig. 6b. 



^In the region oi q> pi the momentum q will destroy the matrix element of high Er.i dijet production 
replacing in (fT2|) the factor 1 /pf by l/q^. 
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Let us finally consider a more realistic situation with Yi = Y2 but p2 < Pi- Recall that 
the true argument of the BFKL amplitude is not rapidity but the momentum fraction x, 
that is actually we have to take Y = ln(l/x). When pi » p2 for the same rapidities 
Yi = Y2 we get in the right ladder the momentum fraction X2 << Xi. In other words, in 
this configuration we may put, in Fig. 5, the recombination vertex just into the cell nearest 
to the left dijet. But then there will be a large \nx (and may be Ing^) interval for the 
evolution of the right ladders (between the dijets on the rhs and the two recombination 
vertices. In other words in Fig.5. we delete only the two 'BFKL blobs' on the Ihs below 
and above the dijet production. Assuming that, in f|T2l) . the total rapidity interval Ytot is 
very large, we may perform first the rapidity integral 

/ '*i'-P[-i'(^(^^) - = ! 2,(,) (29) 

where for the BFKL blobs on the Ihs we have set xif^i) = 0? ^-^^d for /i we put its 
asymptotic value /i = 1/2. Now we close the countor of the /i2 integration around the 
pole x(/^2) — 2x(/i) = 0: this leads to fi2 — 0.18. The same result is obtained for fi^. 
Finally, the integral takes the form 

/pi 
dlng2g2(i-/.2-A.^)^ (30) 

and the major contribution comes from the domain close to upper limit ~ 

A closer look reveals still another detail. In the region of interest, for example in a 14 
TeV pp-coUision at the LHC, we observe in the central region the dijet with pi ~ 20GeV, 
corresponding to x ~ 2pi/y^ ~ 0.003. For such x-values, the anomalous dimension 
observed at HERA is not so large. For x < 0.01 the behaviour of the structure function 
F2(x,g^) can be parametrized as 

F2 = ciq^)x^ (31) 

with c ~ const{q^) and A = 0.048 ± 0.004 [13]. This means that effective anomalous 
dimension /Xe// = A ln(l/x) ~ 0.28 for x = 3 ■ 10~^. This value is still large enough to pro- 
vide the convergence of the integral fl25p in the large q^ domain for the case considered 
above where both recombination vertices are jsut near the dijet production cell. However 
it is not evident that the parametrization fl3Tl) reflects the behaviour of a single ladder. 
At not large q^ the experimentally measured F2 already includes some absorptive effects 
which reduce the growth of F2 with x decreasing and thus leads to a lower value of A in 
comparison with a single ladder contribution. In other words the true value of /Xe// which 
corresponds to a single ladder may be even larger, pushing the characteristic values of q^ 
closer to the (lower) hard scale p\. 



4 Generalizations 

So far we have discussed the effect of two recombinations inside a two-chain contribution: 
one recombination on esch side of the produced jet pairs. Let us first comment on the 
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Figure 7: A recombination of two ladders which allows for diffractive states. 



case where we have no second recombination vertex above the jet pairs: as far as only 
one recombination vertex is concerned, the integration over q is logarithmic. However, 
q runs also through both upper ladders and defines the low momentum scale Qq where 
the evolution starts: a large value of q therefore kills the evolution in the upper ladders, 
whereas a low value prevents the evolution in the lower ladders. Therefore, a single 
recombination vertex is suppressed. 

Next a comment on the color suppression factor (jll). This suppression applies to the 
case when, as illustrated in Fig. 2, there is evolution above and below the recombination 
vertex. As we have discussed before, in a preferred situation we have little or no evolution 
between the recombination vertices and the dijet production vertices. In this case there 
is no need to reconnect, between the two recombination vertices, the four t-channel gluon 
lines to color singlet pairs. As result, the color suppression becomes much weaker.. 

Next we consider the case of more than two chains, say three chains with three pro- 
duced pairs of jets.. In this case a pair of two recombination vertices can be attributed to 
any pair of chains, i.e. we have three possibilities. Similarly, for n chains we have 
possiblities: these counting factors can easily overcome the color suppression factor in (jl]). 
As an example, for = 4, the overall counting factor is already 3/4, and it exceeds unity 
for n>b. 

Finally, we mention another important possibility, related to final states with rapidity 
gaps. Besides the recombination illustrated in Fig. 2 there exists another configuration to 
which our discussion applies. We show this in Fig. 7: Applying our previous discussion, in 
particular the evolution paths illustrated in Fig. 6, we conclude that the momentum scale 
at the upper end of the lower rapidity gap, g^, will be above Qo but not too close to the 
jet momenta p\ = p\. this allows for 'semihard' diffraction and is in qualitative agreement 
with inclusive diffraction seen at HERA. 



5 Conclusions and outlook 

We have studied the possibility of interactions ('recombinations') between two evolution 
chains, which describe double parton scattering corrections in high energy hadron colli- 
sions. We show that, thanks to large anomalous dimensions of the parton distributions in 
the low-x region, such an interaction may occur at small distances within the pertubative 
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domain, provided we consider two recombination vertices (which describe the chain-chain 
interaction) placed relatively close to the 'hard' matrix elements. 

In our double leading log analysis of gluon ladder diagrams we find a competition 
between 'coUinear' logarithms (Inp^) and 'energy' logarithms (ln(l/a;)). Depending on 
the ratio between these logarithms, the major contribution to integrals over the rapidity 
of the two recombination vertices, Y and Y', comes either from the region near the 
protons (Fig.6a) or from the region close to the jet production vertices (Fig.Gb), i.e. these 
integrals have no saddle point somewhere in centre of the available interval. The first case 
(Fig. 6a) corresponds to two independent ladders which do not communicate with each 
other and are described by 'double DGLAP' evolution equations. More interesting is the 
second possibility (Fig. 6b) where the recombination vertices are close to the hard matrix 
elements and thus are entirely in the perturbative region. These configurations may lead 
to nontrivial correlations between the secondaries produced in 'double parton scattering' 
processes. We note that, in Fig. 6b, the rapidities and transverse momenta of the partons 
inside the recombination vertices need not be very close to the 'hard' matrix element: in 
a more or less realistic situation the convergence of the integrals in rapidity {Y and Y') 
and transverse momentm (q) are rather slow, since they are driven by numerically small 
powers oi 1/x and q^. Therefore the particles coming from the recombination vertices 
may still be separated from those produced via the 'hard' subprocess by relatively large 
intervals in rapidity (~ few units) and in the logarithms of transverse momenta. 

The interactions between different ladders discussed in this paper also allow for semi- 
hard diffractive final states (Fig. 7). 

In the case of multiple parton interactions with a larger number of evolution chains 
the suppression of chain-chain interactions caused by the colour factor 1/{N^ — 1) may be 
compensated by the combinatorical factor. For n chains we have n{n — l)/2 possibilities. 
According to naive 'eikonal model' estimates the mean number of chains in proton-proton 
collision at the LHC is < n >~ 5. Therefore the expected probability of multi-chain 
recombinations is not small and may lead both to a noticeable correlation between the 
secondaries in inclusive processes and to 'semihard' diffraction final states. 

All results of this paper are based upon the double logarithmic approximation (with 
fixed as). We consider this as a first step towards a more accurate analysis. Within the 
small-x approach it is possible to go beyond the double logarithmic approximation and 
to reach single logarithmic accuracy (leading In 1/x). Also, a more detailed numerical 
analysis will be needed in order to obtain a more reliable estimate of the importance of 
the recombination corrections addressed in this paper. Both tasks will be topics of future 
work. 

We finally mention an important consequence of our result. In contrast to non- 
interacting multiparton chains which often are modelled within the eikonal approximation, 
corrections due to the recombination of ladder diagrams no longer fit into the eikonal pic- 
ture. This raises the question of the AGK cutting rules which provide a crucial theoretical 
constraint of multiparton corrections. An investigation of this problem is quite important. 
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